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THE  KINETICS  OF  FORMATION  OF  TWO-PHASE  SISTEMS  SEMI  THE 
CRITICAL  POINT 

l 

I 

i 
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/following  is  the  translation  of  an  article  by  la.  B,  j 

Zel'dovich  and  0,  M.  Todes  entitled  "Kiaetika  Obraaovaniya 
Dvukhfaanykh  Sistaa  Vblizi  Kriticheskoy  Tochki  (English  ; 

version  above)  in  Zhurnal  Eksper imental 1 noy  i  Teoret-  | 

icheskoy  Fiaiki  (Journal  of  Experimental  and  S 

j 

Theoretical  Physics)  ;Vol.  10,  No.  12,  194-0,  pp  1441  — 

1445*7  | 

'  » 

The  process  of  stratification  of  a  substance  which  is  in  an 

j 

unstable  state  relative  to  small  changes  near  the  critical  point,  is  i; 

considered.  Of  importance  for  the  stratification  time  is  the  rate  !. 

1 

J 

of  heat  exchange,  and  the  distance  of  heat  transfer  and  the  value  of  j 
the  temperature  gradient  are  estimated  from  the  capillarity  theory.  j 
Limiting  relations  for  the  stratification  times  are  obtained 
by  expanding  the  thermodynamic  quantities  in  powers  about  the  critical 
point. 


V,  T  plane)  the  region  1,  in  which  thermodynamic  equilibrium  ecu- res¬ 
ponds  to  a  single-phase  system,  from  region  II,  in  which  the  system 
breaks  up  in  thermodynamic  equilibrium  into  two  phases  (liquid  and 
vapor). 


Going  over  to  an  arbitrary  point  B  on  this  line,  we  enter  near 
the  latter  into  part  of  region  II,  which  we  shall  call  region  Ila.  j 
In  this  region  the  single— phase  system  is  unstable  relative  to  the 
formation  of  the  new  phase,  which  differs  essentially  in  its  proper-  j 
ties  from  the  initial  phase  Bj  in  Fig.  1,  the  second  phase 
thus  produced  is  represented  by  the  point  B.  The  occurrence  of  an 
essentially  different  new  phase  is  connected,  as  is  well  known,  with  j 

ithe  need  of  consuming  surface  energy}  these  considerations  have  led  j 

Gibbs  to  the  concept  of  the  critical  dimension  of  the  nucleus  of  the 

asw  phase  and  to  an  explanation  of  the  considerable  stability  of  the 

j  ■ ' 

single-phase  system  near  the  boundary  curve  ABODE*  i An  experimental  * 

i 

disturbance  to  super  cooled  vapor  or  to  super  heated  liquid  is  j 

i 

; 

actually  connected  with  such  extraneous  factors  as  the  presence  of  j 

\ 

contaminating  dust  particles,  ions,  etc*  In  the  terminology  of  Gibbs,  j 
Ila  is  the  region  of  states  which  are  stable  under  small  changes  in  j 

Ithe  phase,  but  are  unstable  under  the  occurrence  of  a  new  different  I 


base. 


- K  ibe  geometrical  treatment,  tbs  line~ABG^r7^ seats  a 

projection  on  the  plane  V,  T  of  the  line  of  intersection  of  two 
surfaces  of  chemical  potential.  However,  oa  each  surface  separately,  | 
this  *M"«  of  intersection  does  not  correspond  to  any  singularity,  from| 
which  we  see,  in  particular,  the  stability  of  the  states  which  sre  . 

j 

dose  to  the  boundary  line  ABODE,  relative  to  small  changes  in  tbs  j 

i 

phase.  | 

Near  the  critical  point,  in  which  the  difference  between  the 

existing  phases  in  equilibrium  vanishes,  there  is  a  region  lib  of 
states  which  have  very  remarkable  properties.  Thus,  in  the  region  lib 

|  '  (.w)r^mrt>0* 

!  s 

whereas  on  the  boundary  of  the  region  lib,  along  the  line  FGCIK,  and 
in  particular  at  the  critical  point  C,  we  have 

|  (?;),“«• 

This  region  is  already  the  region  of  states  which  are  unstable  rela¬ 
tive  to  small  changes  in  phase  (Gibbs). 

In  region  lib,  changes  In  the  states,  no  matter  how  small,  will 
increase  ana  will  lead  as  a  result  into  the  separation  of  the  system 

into  two  phases.  At  the  same  time,  in  region  Ha,  small  changes  of 

disappear 

the  state  decrease,  and  completely  as  in  any  stable  state  I.  In 
the  region  Ila  the  stratification  takes  place,  as  already  indicated, 
only  after  the  formation  of  a  nucleus,  which  the  parameters  that 
oharaoteriae  the  state  (specific  entropy,  specifio  volume)  differ  by 


H 


of  finite  volume. 

The  instability  with :  respect  to  smell  changes  in  states, 

i»«idoh  !©,>«•: 

is  dirootly  obvious,  for  the  reduction  in  the  volume  —  compression  - 
leads  to  a  decrease  in  the  pressure  and  a  farther  compression. 

We  consider  below  specifically  the  speed  of  stratification  of 

4  # 

a  system  in  state  Hb,  confining  our  investigation  to  the  only  oast 
of  physical  interest,  thst  of  states  dose  to  the  critical  point. 

.'2,  The  Hiat  Transfer  Determines  the  Speed  of  Stratification 
A  fact  which  is  essential  for  the  question  under  consideration 
is  that  in  the  critical  point,  where  (dp/  dv)i  =  const  reverses  its 
sign  aai  passes  through  aero,  the  value  of  (dp/d  v)g  s  ooaa^  •*“ 
perienoe  no  essential  changes  and  remains  negative. 

Thus,  an  adiabatic  oompreesion  of  certain  elements  of  volume 
and- the  expansion  of  others  will  lead  to  an  inoraaae  in  preesure  in 
t ha, oompre seed  parte  of  the  volume  and  to  a  drop  in  pressure  in  the 
expanding  ones,  and  consequently,  to  a  reverse  equdisation  of  these 
Changes,  relative  to  which  the  single-phase  system  ia  stable  even  in 
region  Hb.  The  instability  indicated  in  the  end  of  the  preceding 
Motion  for  the  states  in  region  Hb  pertains  only  to  isothermal 


proo(tW8«  An  Isothermal  change  in  the  state  is  ooonsoted  with  s 
change  In  the  entropy,  and  consequently,  requires  heat  exchange 
between  the  expanding  and  compressing  elements  of  the  volume. 

The  process  of  stratification  of  a  system  into  two  phases  in 
region  lib  is  thus  connected  with  motion  of  matter  and  with  heat 
transfer.  For  a  single-component  system,  the  first  prooess  does  not 
require  ary  molecular  transport,  and  is  realised  by  molal  motion  of 
masses  and  is  therefore  quite  rapid  oompared  with  the  second  prooess 
—  heat  transfer.  (In  the  case  of  the  critical  point  for  binary 
mixtures  (far  example,  water-phenol),  the  appearance  of  which  is 
connected  with  the  separation  of  a  homogeneous  solution  into  two 
phases  of  different  composition,  one  can  expect  that  the  limiting 
factor  will  be,  to  the  contrary,  not  heat  transfer,  but  the  prooess 
of  transfer  of  matter,  realised  by  molecular  diffusion.) 

I  Therefore,  in  the  process  of  disintegration  of  the  states  in  region 
lib,  one  can  consider,  the  pressure  to  be  instantaneously  equalised 
in  practice.  Further  calculations  of  the  heat  transfer  will  involve 
the  specific  heat  at  oonstant  pressure.  The  instability  of  the  system 
relative  to  small  changes  in  state  is  reflected  in  the  singular 
behavior  of  the  apeoifio  heat.  According  to  the  relation 

the  ohange  in  the  sign  of  (df>/&v)T,  passing  through  sero,  implies 
a  change  in  the  aign  of  Cp  with  transition  through  a  CO  .  It  is 


•••A  therefore  that  In  ths  Region  IBs  tb#  tptoiflo  teat  Cp  is  nega¬ 
tive.  The  change  in  the  sign  of  the  speoifie  haat  in  the  aquation  of 
i  haat  conduction  ie  formally  equivalent  to  a  change  in  the  sign  of  the 
time;  instead  of  ordinary  exponentially  decreasing  solutions,  des¬ 
cribing  the  equalisation  of  the  temperature ,  vie  obtain  exponentially 
Increasing  temperature  differences,  whioh  lead  to  a  stratification 
into  two  phases  with  subsequent  equalization  of  temperature  after 
the  limiting  line  FGCIK  if  reaofaed. 

*  The  amplitude  of  the  primary  fluctuation,  the  exponential 
growth  of  which  lead#  to  stratification  into  two  phases,  will  enter 
into  the  expression  far  the  time  of  stratification  logarithmically, 
and  therefore  is  of  no  importance  in  ths,  approximation  adopted  in  qur 
paper. 

Let  os  expand  the  free  energy  7,  the  pressure  p,  and  the  thermo- j 
dynamic  potential  0  in  powers  of  the  distance  to  the  orltioal  point  j 
(in  variables  V  and  T)  about  the  critical  point.  j 


Vk)~a(T -  TMV~  Vk)+ 
-*-b(T~Tk)(V—  Vkf  +  d(V-  V? 

Vt)[2b(T-  Tk)-t-4d(V —  Vf\ 
<b=:FQ-^-pVk — 6(7’ —  Tk)(V-  VkY-3d(V~  VJ\ 


(2) 

P) 

<4> 


<*’  i 


In  the  7,  7  plane,  the  equation  of  tha  boundary  line  ABCDB  is 
leter mined  from  the  conditions  Pi  s  P2  ao^  01  s  02*  yi«ld® 


=vk±yi 


H&rZi. 

2 d 


(5) 


7 


An  estimate  of  the  time  of  stratification,  from  a  consideration 
of  the  equation  of  heat  conduction 


oC^-KAT 


f 

(  ^  a  1/V  is  the  density  and  K  is  the  coefficient  of  heat  conduction)  ,i 

i 

requires  the  determination  of  the  characteristic  dimension  /,  i.e.,  ! 

the  geometric  dimensions  of  those  regions,  on  which  the  homogeneous  j 
initial  system  is  stratified.  The  sought  time  is,  in  order  of  magnl-  j 


Tfaa  stratification  is  the  £ K  the  farther  away  we  are  from  the 
oritioal  point,  where  °°  j  it  ia  also  the  faster,  the 

smaller  the  critical  dimension  /.  An  estimate  of  the  minimum  dimen¬ 
sion  Jf  requires  the  has  Of  additional  considerations. 

We  shall  use  the  classical  theory  of  capillarity  near  the 
oritioal  point.  The  law  of  variation  of  surface  tension  on  the 
boundary  between  liquid  and  vapor  ia  near  the  oritioal  point 

a=a(n-7f*  (10) 

is  derived  from  this  theory  from  ideas  oonoerning  the  term  in  the 

[ 

thermodynamic  potential,  dependent  on  the  density  gradient.  Erom 
symmetry  considerations  it  follows  that 

<P=$(r,  »)-*-ns<V«)s.  (ll) 

It  is  obvious  that  the  added  term  insures  a  continuous  ohange-  in 
state  on  the  separation  boundary  of  the  two  phases.  The  effect  of 
width  of  the  separation  boundary  is  determined  from  the  condition 
of  minimum  of  thermodynamia  potential  of  the  system  Y  -  J  $  dm 

O 

(integral  over  the  entire  mass  of  the  substance  M). 

Everywhere,  except  on  the  separation  boundary  between  the 

two  phases,  V  *  or  V  *  V2  fees  (5fl  J  0  *  0(V>  *  ^V2,T)  =  00- 
In  a  surface  layer  having  a  volume  /  there  is  contained  a  mass 

.41  .  *  *  «  .  * 

SljL/V kj  we  put  there  $(V,  T)—$ >  T)==*<-v*  * 

:  [cp.  (4),  (5)1 

a  *  -  — 

Denoting  f(js  0qM,  we  obtain 


(12) 


vhere  -O.  is  the  total  area  of  the.  separation  boundary  of  the  two 
phases.  From  this  we  obtain 


t^rtfzrry  0 ~ T?~  a(Tk  —  Tf 


By  comparison  of  the  obtained  limiting  laws  (5),  (7)  and  (10)  with 
the  experimental  data  on  surface  tension,  specific  volume,  and 

%  %  *  *  «  t 

specific  heat  near  the  critical  poiatj  we  obtain  the  coefficients  n> 
b,  and  d  which  are  involved  in  the  calculatione  The  relaxation  time 
of  interest  to  us  is 


1  n*  a 1  7* 

%  ~ &r*  r(ft  -  f)2i(n  -  v,)- ' 


Assume  that  we  obtain  a  state  in  the  region  lib  by  adiabatic 
expansion  of  the  gas  such  that  the  adiabat  passes  exactly  through 
the  oritical  point.  Hare  we  have  hear  the  critical  point  along  the 


adiabat 


r»  —  vt~a(Tr-  T),  1 

d(V—  Vf<r-2b{Tt~  T)i 


and  the  stratification 


time  of  the  system  is 


~(7>-  T)1. 


4.  Estimate  of  Order  of  Magnitude  of  the ^Stratification  Time 
Accurate  to  numerical  factors,  let  ua  express  the  coefficients 
contained  in  (14)  in  terms  of  the  critical  constants  of  the  substance. 
as  can  be  seen  from  the  dimensionalities 


/o 


j  *  •'TO* 


Furthermore,  we  pot  approximately 


#s^ 


uhere  €  q  is  the  surface  tension  with  T  ^  Tjc»  Substituting  (17) 
and  (18)  in  (14),  we  g«t  - 


Far  water  in  absolute  units  &o  ■=£  80,  Tk 
and  K..S  10^,  we  obtain 


_ 2-10-9 

r=f?T=Tp 


S'  60CP  K,  pn  «  2  x  10* 


It  nn#  be  noted  here  that  in  the  derivation  of  (20)  we  omitted 
everywhere  the  numerical  faotors  whioh,  upon  accumulation,  can  con¬ 
siderably  inorease  .  If  we  use  in  the  determination  of  a,  b,  and 
d  the;  expansion  of  the  Van  der  Waala  equation  near  the  critical 
point,  we  obtain  far  formulas  (19)  and  (20)  a  ,  *  factor 

1216  x3'3  &  2,000. 

A  further  refinement  of  the  width  /  and  determination  of  the 
time  V'  can  yield  additional  numerical  faotors. 


H 


5.  Unattainable  States 

¥e  hot®  that  far  inside  the  region  Hb  one  can  also  imagine 
states  for  which  not  only  (^P/*V)T  0,  but  also  (3?/*V)s  >  O.j 

However,  such  states  are  not  only  unstable,  but  unattainable  j 

at  all.  Actually,  according  to  (1),  on  the  very  boundary  of  this 
region  the  speolfio  heat  Cp  vanishes,  and  aooording  to  (9)  the  j 

stratification  time  also  vanishes.  Thus,  as  one  penetrates  deeper  j 

into  the  region  II,  the  time  of  stratification  beoomes  shorter  and  j 
shorter  and  on  the  line  <$IS  the  stratification  should  occur  prac-  .  j 

,  *  i 

tically  instantaneously.  j 

In  the  case  when  the  single-phase  system  is  accurately  das-  :  j 

■  '  .  '  1 
cribed  by  the  Van  der  Waals  equation  of  state,  the  position  of  the  j 

curve  QRS  was  calculated  in  1937,  at  our  suggestion,  by  student  ! 

I 

A.  S.Viglln  of  the  Leningrad  Industrial  Institute.  It  was  found  that  ; 
the  minimum  R  of  the  curve  QRS  is  located  precisely  under  the 

1 

maximum  of  the  boundary  curve,  l.e., 

i 

"  '  •  *  I 

:  V,L=tV«  r*~T*T&Z (21)  j 
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